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Abstract 

The degrading effect of RF impairments on the performance of wireless communication systems is more 
pronounced in MIMO-OFDM transmission. Two of the most common impairments that significantly limit the 
performance of MIMO-OFDM transceivers are IQ-imbalance and phase noise. Low-complexity estimation and 
compensation techniques that can jointly remove the effect of these impairments are highly desirable. In this paper, 
we propose a simple joint estimation and compensation technique to estimate channel, phase noise and IQ-imbalance 
parameters in MIMO-OFDM systems under multipath slow fading channels. A subcarrier multiplexed preamble 
structure to estimate the channel and impairment parameters with minimum overhead is introduced and used in 
f\j 1 the estimation of IQ-imbalance parameters as well as the initial estimation of effective channel matrix including 

^ . common phase error (CPE). We then use a novel tracking method based on the second order statistics of the inter- 

^ carrier interference (ICI) and noise to update the effective channel matrix throughout an OFDM frame. Simulation 

. results for a variety of scenarios show that the proposed low-complexity estimation and compensation technique can 

efficiently improve the performance of MIMO-OFDM systems in terms of bit-error-rate (BER). 

I. Introduction 

c/j \ The reliability and throughput offered by orthogonal frequency division multiplexing (OFDM) and multi-input- 
^ ; multi-output (MIMO) techniques has made MIMO-OFDM a clear choice for future generation wireless communica- 
tion systems [fl]]. However, MIMO-OFDM systems are highly sensitive to various radio frequency (RF) impairments 
due to the imperfections in analog front end components. Among various RF-impairments, IQ-imbalance and phase 
^ ■ noise are considered more detrimental than others and can cause severe performance degradation in MIMO-OFDM 
\Q ■ systems. 

Due to flexibility and low complexity, direct conversion and low intermediate frequency (IF) receivers are preferred 
to superheterodyne receivers. These receivers provide perfect image attenuation eliminating the need for image reject 
filters J2). In practice, however, due to manufacturing imperfections, there is always phase and amplitude mismatch 
^ \ in the inphase (I) and quadrature (Q) branches of the IQ-receiver commonly known as IQ-imbalance 0- The 
compensation of IQ-imbalance in MIMO-OFDM has gained a lot of attention recently HI. Il5l. 

Phase noise is another major impairment caused by non-ideal oscillator which introduces random phase fluctu- 
ations at the oscillator output resulting in distortion of the signals Q. The effect of phase noise is more severe 
in OFDM systems and is studied extensively in the literature |6), Q. These works state that the effect of phase 
noise can be broken down into common phase error (CPE) and inter-carrier interference (ICI) with CPE being more 
prominent for systems with small phase noise. Performance analysis and compensation techniques are proposed for 
OFDM and MIMO-OFDM systems in the literature and can be found in Q-(9l. 

Joint IQ-imbalance and phase noise compensation techniques for single-input-single-output (SISO) systems have 
been previously discussed in iflOl and ifTTl . However, joint treatment of these impairments for MIMO-OFDM 
systems have not been studied extensively. To our knowledge, only recently in Ifl2l . a joint estimation and 
compensation scheme for MIMO-OFDM systems assuming quasi-static channel is proposed. The authors first 
develop a compensation scheme for SISO systems and generalize it to the MIMO case. The proposed method 
requires more training symbols with increasing number of transmit-receive antennas and the complexity of the 
scheme is high due to the required maximum likelihood (ML) estimation. In this paper, we first introduce a 
subcarrier multiplexed preamble structure to estimate the channel and impairment parameters with small overhead. 
The IQ-imbalance parameters as well as the initial estimation the effective channel matrix including common phase 
error (CPE) are then estimated on the preamble. Using scattered pilot symbols and second order statistics of ICI 
and noise, we then propose a novel tracking method based on the minimum mean square error (MMSE) to update 
the effective channel matrix in the rest of the OFDM frame. Simulation results for different phase noise and IQ- 
imbalance parameters show that the proposed low-complexity estimation and compensation technique can efficiently 
improve the performance of MIMO-OFDM systems in terms of bit-error-rate (BER). 
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Throughout the paper, upper-case bold letters are used for matrices, lower-case bold letters for vectors and 
scalers are represented by italic letters. Operators *, and indicate convolution, element-wise multiplication 
and element-wise division respectively. {.}*, {.} T , {-} H , {•} and {.}* represent conjugate, transpose, hermitian 
transpose, pseudo-inverse and conjugate minor operation, respectively. Ijv, Oat represent N x N identity and 
zero matrix. E{.}, arg{.}, abs{.} represent expectation, argument and absolute value of the elements inside the 
parentheses, respectively. 

II. System Model 

A. Signal Model 

The block diagrams of typical MIMO-OFDM transmitter and receiver with M t and M r antennas are shown in 
Figs. 1 and 2, respectively. At transmitter, the data stream is first converted to parallel streams and passed through 
iV-point IFFT block after adding pilot symbols. The resulting parallel streams are converted back to serial and a 
cyclic prefix (CP) of length N cp (L < N cp ; where L is the number of resolvable channel paths) is then added to 
the signal stream. The resulting signal is up-converted with suitable modulation and transmitted. Let s p m be the 
m th N x 1 frequency domain OFDM symbol to be transmitted from the p th transmit antenna whose time domain 
representation is given by s p m = T~ (s P)77l ) where J 7 " 1 is the inverse Fourier transform operation. The vector 
that contains all the N + N cp + L — 1 discrete time domain received signal samples at the q th receive antenna is 
then given by 

M t 

= h<j,p * §p,m + w <? (1) 
p=l 

where h q:P = (h qtP (0), hq tP (l), ...,h q>p (L — 1)) T is an L x 1 vector representing the channel impulse response of 
the multipath channel from the p th transmit to the q th receive antenna for the m th OFDM symbol, s p m is the 
transmitted discrete time signal vector from the p th transmit antenna and w q is an (N + N cp + L — 1) X 1 vector 
containing the samples of complex circularly symmetric additive white Gaussian noise (AWGN) at the q th receive 
antenna. At each receiver branch, the received signal is down-converted, sampled and have the cyclic prefix removed 
(see Fig. 2). The signal samples are then converted into parallel streams and fed into FFT block with size N. The 
equivalent frequency domain received signal model at the k th subcarrier (k = 0, 1, .., N — 1) after FFT block can 
be written as 

r m (k) = H m (k)s m (k) + rj m (k) (2) 

where r m (k) is the M r x 1 received signal vector, H m (fe) is the M r x M t channel frequency response matrix, 
s m (k) is the M t x 1 transmitted signal vector at the k th subcarrier of m th OFDM Symbol. r) m {k) is the M r x 1 
frequency domain noise assumed to be complex zero mean AWGN with variance cr^, i.e., r) m (k) € C7V(0, cr^). 
Each (q,p) element of the channel frequency response matrix is given by Hm P (k) = Yln=o hq,p( n ) e ~ J= ^\ 

B. IQ-Imbalance Model 

The time domain oscillator output with frequency independent IQ-imbalance (see Fig. 2) can be written as O, 

M 

c(t) = cos (u) c t) + je sin(ujct + 9) 

= K ie jujJ + K 2 e~ jUct (3) 

where oj c , 6 and e are the carrier frequency, phase mismatch and amplitude mismatch at the receiver, respectively, 
and K\ = (1 + ee~^ e )/2, K% = (1 — ee^ e )/2. It is assumed that each receiver branch uses separate oscillators and 
IQ-imbalance parameters are different for each receiver branch. The received signal in frequency domain after FFT 
block can then be written as 

x m (fc) = Kir m (fc) + K 2 r#(A;) + Vm {k) (4) 

where x m (k) and rj m (k) are the M r x 1 received signal and noise vectors at k th subcarrier of the m th OFDM 
symbol, respectively, and # operation represents conjugate mirror operation i.e., r|(fc) = r* m (—k). The matrices 
Ki and K 2 are given as Ki = (I + ee~ j0 )/2 and K 2 = (I — ee +j0 )/2 where e = {e-J*^ and 6 = {(9*}^ are 
the amplitude and phase mismatches associated with each of the receiver branches. 
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Fig. 1. MIMO-OFDM transmitter with M t transmit branches 
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Fig. 2. MIMO-OFDM receiver with M r receive branches affected by IQ-imbalance and phase noise 



C. Phase Noise Model 

In the presence of phase noise, the local oscillator output is given by c(t) = e J '( We * + ^w) where 4>{t) is the random 
phase fluctuations of the local oscillator. Phase noise process modeling has been extensively studied in the literature 
for both phase locked loop (PLL) based and free -running oscillators Q, Ifl5l . Here, we consider a free-running 
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oscillator model. The phase noise for this oscillator can be modeled as Wiener process with independent Gaussian 
increments such that 4>{n + 1) = 4>(n) +u(n), where <f>{n) is the discrete Brownian phase noise process at n th time 
instant, u(n) G A/"(0, 47r/3T s ) with j3 being the 3-dB linewidth of the phase noise spectrum and T s the sampling 
time period of an OFDM symbol. The time domain received signal corrupted by phase noise at the q th receiver 
can then be written as 



e j<P(0) g#(l) _ e j<f(N+N cp +L-l) 







M t 

^2 h q,P * + W <? (5) 

p=l 

where e J< ^ n ) represents one realization of phase noise process for the q th receiver branch. Equivalently, the received 
signal in frequency domain can be written as 

N-l 

r m(k) = ^2 ~~ k)H m (i)s m (i) + r) m (k) (6) 

i=0 

where & m (i) is the M r x M r diagonal matrix of phase noise frequency domain samples for the i th subcarrier 
(i = 0, .., N — 1) of the m th OFDM symbol. The effect of phase noise for the received symbol can be separated 
into two terms, the CPE (for i = k) and ICI (for i ^ k) terms, as Q 

r m (k) = ® m (0)H m (k)s m (k) + C(k) + Vm{k) (7) 

CPE ICI 

The first term in |7J), m (O), representing CPE is an M r x M r diagonal matrix showing the rotation of the received 
symbols by a common phase while the second term, C(^)> i s the ICI term caused due to phase noise components 
from all other subcarriers at k th subcarrier given by C(^) = Yli^oi^k ®m{i — k)H m (i)s m (i) . CPE is common 
for all subcarriers within one OFDM symbol but varies from one OFDM symbol to the other. Studies show that in 
practice CPE is more harmful to the system performance than ICI at realistic phase noise regimes where ^ C 1 /NT S 
0. 



D. Combined Phase Noise and IQ-imbalance Model 

The combined effect of the phase noise and IQ-imbalance on the received symbol in frequency domain can be 
modeled by combining (01) and {/J as 

x m (fc) = K 1 m (O)H m (fc)s m (fc) 

+ K 2 ®* m (0)H m #(k) Sm #(k) 

+ <(k)+V m (k) (8) 

(H]) reveals that the desired signal at particular subcarrier is affected by the joint effect of multiplicative and additive 
impairments due to IQ-imbalance and phase noise from all the subcarriers. We will show, in the results section, that 
a joint compensation of these impairments should be considered in order to achieve acceptable system performance. 

III. Initial Channel and Impairment Parameters Estimation 
A. OFDM Frame Format 

In typical MIMO-OFDM systems, such as those based on IEEE 802.1 la/n standard, the transmitted OFDM 
symbols are divided into training and data symbols lfl6l . Training OFDM symbols are further classified into short 
training symbols and long training symbols also called as preamble. Furthermore, data type OFDM symbols are 
usually loaded with pilot symbols for fine tuning the estimation of channel and impairments. The null subcarriers 
in the short OFDM training symbol can also be used to estimate the noise plus ICI correlation matrix as given in 
j9ll . For L OFDM short symbols with N nu u null subcarriers, this correlation matrix \l/ can be estimated as (9[ 

L N nuU 

t=l k=l 
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B. Effective Channel and IQ-Imbalance Estimation 

The accuracy in estimating channel and impairment parameters at the preamble stage of OFDM transmission 
largely depends on the OFDM training (preamble) symbols design. Due to the presence of IQ-imbalance, interference 
from mirror subcarriers should also be considered in preamble design. On the other hand, utilizing more OFDM 
symbols for channel estimation at the preamble stage results in poor estimation performance due to different CPE 
at each training symbol ifTTl . In this paper, we use a modified subcarrier multiplexed preamble structure |[T8l with 
minimum overhead that is more robust against the varying nature of CPE. We define matrices Ti and T2 that 
contain M t training symbol vectors of size N x 1 to be transmitted from M t transmit antennas as 

Ti = {7 £i,7 £2,-, 7 °£mJ (10) 

T2 = {7°£i,7°4-,7°£mJ (H) 

where 7 is a N x 1 vector of known training symbol that can be optimized according to the system requirements 
and the N x 1 vectors £ p and £ p corresponding to p th transmit antenna are given as 

n p _i{[i, o Mt -i, 1, o Mt -i, 1, o Mt -i]} T (12) 

-i p {k) for ke {1,2,., f -1} 

(13) 

£ P (fc) for ke {-f,-f + 1,..,-1} 

In (fT2l . n p _i represents the cyclic shift operation over p — 1 samples and 0m 4 -i denotes an all-zero vector. Fig. 
3 shows the preamble structure for the special case of 2 transmit antennas. With the above proposed preamble 
structure and without considering ICI and noise, the received signal vector for the two consecutive training symbols 
can be written as 

tpi(k) = K 1 H pre (k)t 1 (k) + K 2 H* e {k)4(k) (14) 
il> 2 (k) = Ki(H pre (fc) + S)t 2 (k) 

- K 2 (H# e (fc) + S*)tf (k) (15) 

where t\(k) and t 2 (k) are the k row of Ti and T2 with each of them containing only one nonzero element Xi(k) 
and X 2 (k) respectively. H pre (A;) = pre (O)H(/c) and H.f re (k) = 0* re (O)H#(A;) represent the effective channel 
matrix at the k th and —k th subcarriers with pre (O) being the CPE associated with the training symbols. 6 is the 
M r x 1 vector representing error term caused due to different CPEs affecting two consecutive training symbols. 
This error term is the same for all the subcarriers. We define two vectors, Xa{k) and Xb{k) as 

Xa(k) = {Mk) + ^(fc)}/2Ai(fc) = Kifi„(A0 + p (16) 
Xb (k) = Wi(k) - Mk)}/2X*(k) = K 2 hf(k) - p (17) 

where h p (k) = G pT . e (0)h p (A:) and hf(k) = @* re (0)hf(k) represent the p th column vector of H pre (A;) and 
HpVe(fe), respectively and p = (Ki<5 - K 2 S*)/2. Using £[6]> and ([17]) and the fact that K 2 = Im t ~ K^, we can 
approximate h p (k) as 

e(fc) = Xaik) + X f{k) = hp(fe) + p - p* (18) 

From (fT8T ). the channel vector corresponding to the k th subcarrier (k = kMt +p, k = 0, 1, .., [N/Mt]) of the p th 
receiver branch is estimated with very small error, p — p* = | (6 — S*), that does not depend on the IQ-imbalance 
parameters and is caused solely due to CPE difference between two consecutive long training symbols. It is shown 
later in section IV that this error can be compensated for by applying a tracking algorithm. The channel vectors 
for the rest of the subcarriers can be estimated using interpolation or transform domain iterative techniques |[T9l , 
||2U1 . To estimate the IQ-imbalance parameters, for k = 0, N — 2, we define 

a(k) = e(k)-e(k + l) 

= h p (k)-h p+1 (k + l) (19) 

P(k) = Xa(k) - Xa(k + 1) 

= Ki^AO-Kihp+itfc + l) (20) 
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(b) Second OFDM Training Symbol 
Fig. 3. OFDM preamble structure for 2 transmit antennas 



The phase and amplitude imbalance parameters are then calculated as 

6{k) = -arg{2(/3(fc)0a(fc))-l Mr } 
e(k) = abs{2(/3(k) ct(k)) - l Mr } 



(21) 
(22) 



where is the M r x 1 all one vector and represents element-wise division operation. The obtained imbalance 
parameters can be averaged over all subcarriers to get a better estimate. To further increase the accuracy of the 
estimates and because of the slow varying nature of IQ-imbalance parameters, these parameters can be averaged 
over several OFDM symbols (or frames). 



IV. Effective Channel Tracking and Data Detection 

We assume quasi-static channel such that the channel response can be considered fixed within an OFDM frame. 
The CPE term, however, changes from one OFDM symbol to the other and needs to be estimated and compensated 
for each OFDM symbol. Assuming there are r pilot subcarriers in m th OFDM symbol from j transmit branch 
with pilot values d m j(l), d m j(2), d m j(r), it is possible to update the effective channel matrix estimates in 
(fT8l) at data transmission stage as follows. The M r x 1 received pilot signals at I th and — I th subcarriers of the m th 
OFDM symbol where I G 1, 2, ..,r can be written as (without considering ICI and noise) 

x m (Z) = K!T m y m (Z) + K 2 T^y#(/) (23) 
x# (I) = KtT* m y*(l) + K* 2 T m y m (l) (24) 

where, T m = diag{T mi i, .., T m ^M r } represents the CPE updating matrix such that T m ^ q = Q m ,q{0) / ®pre,q{0) and 
Ki = diag{k^\ .., k^ ) and K 2 = diag{k^\ .., k^ ) are the IQ-imbalance matrices. We also define vector y m (/) 
as y m (0 = {y m ,i(0i ••) ym,M r (l)} T such that y m (Z) = H pre (/)d m (/) where fl pre (l) is the estimated channel matrix 
using training symbols as given in section III and d rn (l) = {d m ^i{l), .., d m: M t (0} T i s tne vector of transmitted 
pilots from different transmitter branch at I th subcarrier. 

Algorithm Q] summarizes the procedure to compute the updating parameter matrix to compensate for varying CPE 
for each OFDM symbol. As the CPE is the same for all the subcarriers, it is possible to get a better estimate of 
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Algorithm 1 Calculation of updating parameter matrix Y 

Input -> Ki,K 2 ,H pre (fc),*,r,x m (Z),x#(Z) 
Output — > Y m 

• Initialize — > q — 1 
while q < Mr do 

Fug l\.q 
L.*(2) l,*(2) 

* Initialize — ¥ I = 1 
while Z < r do 

1. Set z q (l) = [a;m, 9 (0,^m,g(0] T 

2. Compute y m , q {l) and y£,,(Z) 

3. Compute C,(Z) = y m , q (Z)Xi i9 +y^, g (Z)X 2 ,, 

4. Compute Q 9 (Z) = (Cf (Z)C,(Z) + * 9 , g I 2 )~ 1 Cf (Z) 

5. Compute <^ q (Z) = Q q (l)z q (l) 

6. Z = Z + 1 

- Compute tp q = 1 £[ =1 £,(Z) 

- Compute T m , 9 = <£ 9 (1) + j(p q {2) 

- q=q+^ 

• Set T m = diap-fTm,,}^ 



Set Xi„ 



jL*(2) 7, 



(i) 
<? 

(2) 



X 
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the updating matrix by averaging it out over all pilot subcarriers. After finding the updating matrix, the transmitted 
data vector at the k th subcarrier, s m (k), can be estimated using H m (/c) 

2 



T m H pre (A;), Ki, K2 and solving 



min||s m (fc) - A(k)x m (k)\\' (25) 

where s m (k) = {s m (k), and x m (fe) = {x m (i),xf (fc)} T are the estimated and received signal vectors, 

respectively. For zero forcing (ZF) and minimum mean squared error (MMSE) based receivers 

A ZF (k) = (W H (k)W(k))- 1 W H (k) (26) 
A M MSE(k) = (W H (k)W(k)+Ky 1 W(k) H (27) 



where 



W(jfe) 



KiH m (/c) K 2 H#(A;) 



K* 2 U m {k) KlH^ik) 



(28) 



and R = * ® I Mr . 

The complexity of our proposed algorithm can be analyzed in terms of the required complex additions and 
multiplications. The channel estimation in (fT6l)-(fT8T) only requires 3N complex additions and 2N complex mul- 
tiplications. Also, the IQ-imbalance parameters estimation in (fT9l)- (l2"2l requires 2A^ complex additions and 2N 
complex multiplications. The estimation process does not involve any matrix inversion. The tracking algorithm 
involves inversion of a 2 x 2 matrix for rM r times. It also involves extra 13rM r complex multiplications and 4rM r 
complex additions. The scheme proposed in |[T2l is more complex than our proposed scheme since it involves 
higher number of complex multiplications. Note that in 11121 . the analysis is done only for SISO systems while the 
computational complexity increases exponentially for by the number of transmit/receive antennas. 



V. Simulation Results 

To validate the proposed scheme, a MIMO-OFDM system with spatial multiplexing transmission scheme |[T3l 
using 16-QAM modulation format is considered. The channel is assumed to be 7-tap multipath Rayleigh fading with 
exponential power delay profile for each transmit-receive link. We consider the 802.11a/n standard frame format 
in which each OFDM symbol is composed of 48 data, 4 pilot and 12 null subcarriers |[T6l . Cyclic prefix longer 
than the channel length is added at the beginning of each OFDM frame. The OFDM sampling time is considered 
to be T s = 0.05/js. One short training symbol and two long training symbols are transmitted according to (fTOl ) 
and (flTT ) in the beginning of each frame. One OFDM frame consists of 50 OFDM symbols (including training and 
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Fig. 4. MSE of channel estimation for 2 x 2 and 4x4 MIMO channels with IQ-imbalance of (5°, 10%) and different phase noise linewidths 
at each receiver branch 



data symbols) during which the channel is considered to be fixed. The IQ-imbalance parameters are chosen to be 
5°, 10% for each receiver branch. 

Fig. [4] shows the mean-squared error (MSE) of effective channel estimation (CE) for 2 x 2 and 4x4 MIMO 
channels with different phase noise linewidths, j3. For channel estimation utilizing interpolation, spline interpolation 
technique is used and for iteration technique, the algorithm presented in (I9\i is used with 50 iterations. As shown, 
an increase in the linewidth of the phase noise increases MSE for both techniques as it introduces higher ICI. Also, 
there is little gain in terms of MSE in using iteration technique for 2 x 2 MIMO system compared to interpolation 
technique; however, this gain is larger in the case of 4 x 4 MIMO. The flooring effect in MSE at high SNRs is 
because of the error propagation caused due to the initial channel estimation error. 

In Fig. [5j the MSE performance of the estimation of IQ-imbalance matrix (Ki) using the metric E{||Ki — Ki|| 2 } 
is shown for a 2 x 2 MIMO channel with different phase noise linewidths. The estimated IQ-imbalance matrix, 
Ki, is calculated by averaging the phase and amplitude imbalances over all the subcarriers at the preamble stage 
and also over two consecutive OFDM frames. As shown, the IQ-imbalance parameters estimation performs very 
well even at high phase noise linewidths. 

Figs. [6] and [7] show the performance in terms of bit error rate (BER) with phase noise linewidth of /3 = 5 
KHz and IQ-imbalance of 5°, 10% at each receiver branch. As shown, the BER curves of the uncompensated and 
partially compensated systems have error floors at mid-range SNRs due to the combined effect of phase noise and 
IQ-imbalance. The performance of the proposed scheme, however, shows significant improvement compared to the 
ones that only compensate for IQ-imbalance or phase noise. 

VI. Conclusion 

In this paper, a joint estimation and compensation scheme to mitigate IQ-imbalance and phase noise in MIMO- 
OFDM systems was presented. Assuming a suitable preamble structure, the channel and IQ-imbalance are initially 
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Fig. 5. MSE of IQ-imbalance matrix (Ki) estimation for a 2 x 2 MIMO system with IQ-imbalance of (5°, 10%) and different phase noise 
linewidths at each receiver branch 



estimated. Then, using pilot symbols in data stage of OFDM transmission, a method to update the effective 
channel estimates including CPE was discussed. The effective channel estimates can then be used to detect the 
transmitted data using detection schemes such as ZF and MMSE. The performance of the proposed scheme for 
spatial multiplexing in MIMO-OFDM systems was verified through simulation and it was shown that it provides 
significant performance improvement in terms of BER with reasonably low computational complexity. 
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Fig. 7. BER performance of an uncoded 16-QAM, 4x4 MIMO-OFDM system with /3 
receiver branch 
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